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Abstract 

We discuss the so-called Schrodinger problem of deducing the microscopic 
(basically stochastic) evolution that is consistent with given positive bound- 
ary probability densities for a process covering a finite fixed time interval. 
The sought for dynamics may preserve the probability measure or induce its 
evolution, and is known to be uniquely reproducible, if the Markov property 
is required. Feynman-Kac type kernels are the principal ingredients of the 
solution and determine the transition probability density of the correspond- 
ing stochastic process. The result applies to a large variety of nonequilibrium 
statistical physics and quantum situations. 



1 Feynman-Kac kernels and time adjoint pairs of 
parabolic equations in the description of ran- 
dom dynamics 

The Schrodinger problem Q of reconstructing the "most likely" interpolating dy- 
namics which is compatible with the prescribed input-output statistics data (an- 
alyzed in terms of nowhere vanishing boundary probability densities) for a pro- 
cess with the time of duration T > 0, can be given a unique solution, 0. For 
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this purpose, it is necessary to define a suitable transition probability m(A, B) = 
J A dx J B dym(x,y), mapping among Borel sets A — > B in time T, so that: (a) the 
bi-variate density m(x,y) has the boundary data Po(A), Pt(B) as its marginals for 
all A and B, (b) m(x,y) has the product form with a certain strictly positive and 
jointly continuous in all variables kernel as a factor, |2|, §. If the respective kernel is 
associated with a strongly continuous dynamical semigroup, then the interpolating 
process is Markovian, || 

The major issue, always to be addressed is: to specify under what circumstances 
(possibly phenomenological, like in case of the boundary density data) the kernel 
can be selected as appropriate, in reference to a concrete physical situation. Clearly, 
the obvious and natural candidates with a direct physical appeal are the familiar 
Feynman-Kac kernels, || [|, [5|, ^j. 

In the physical literature a standard arena for the usage of Feynman-Kac ker- 
nels, and the related Feynman-Kac representation formula for solutions of parabolic 
partial differential equations, is either the Euclidean quantum theory [IT], |3|, [12] , 
or the statistical physics of nonequilibrium phenomena. For example, the Fokker- 
Planck equation, with its non-Hermitian Markov generator, is casually mapped into 
the parabolic evolution problem, whose (semigroup) generator is selfadjoint. In- 
directly 0, [15|, ^ 0, through the Cameron-Martin formula, the Feynman-Kac 
kernels appear as an important tool of the so called stochastic analysis of measures 
and related stochastic processes. It is not accidental, since probability measures 
and their densities are involved in each of the considered frameworks and studying 
dynamics in terms of densities |jl6j [L7j is a theory with much broader, both physical 
and mathematical range of applications, than indicated above. 

The Schrodinger equation and the generalized heat equation, which is basic for 
the original |18| Kac formula derivation, are connected by analytic continuation in 
time. (We shall proceed in the notation appropriate to problems in space dimension 
one, although the main body of our arguments is space dimension independent). 
For V = V(x),x G R, bounded from below, the generator H = —2mD 2 A + V is 
essentially selfadjoint on a dense subset of L 2 , and the quantum unitary dynamics 
exp(—iHt /2mD) is a final result of the analytic continuation procedure for the holo- 
morphic |0| semigroup exp(—Hcr/2mD), a = s + it, s > 0, s — > 0. Here, by equating 
D = h/2m, the traditional notation is restored. Since the unit ball in L 2 is left in- 
variant by the unitary dynamics, and the Born statistical interpretation postulate 
assigns to each normalized function ip(x,t) = [exp(—iHt/2mD)ip](x,0) a probabil- 
ity measure fi(A) = J A p(x,t)dx < 1 with the density p(x,t) = ip(x,t)ip(x,t), we are 
quite naturaly facing the problem of the existence of the random dynamics (stochas- 
tic process) which is compatible with the given time evolution of p(x, t), or preserves 
the measure in the stationary case. 
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Let us emphasize that it is notour goal to propose any probabilistic "derivation" 
1 21] of quantum theory. Rather, we take seriously the Born postulate and attempt 
to draw consequences of this assumption. Its impact is not merely conceptual: the 
mathematical structure of the theory is affected by submitting the quantum unitary 
dynamics to the methods of stochastic analysis, appropriate for any standard prob- 
abilistic problem. Quite irrespectively of whether we deal with essentially classical 
or quantum phenomena, and whether they are intrinsically random or have merely 
a random appearence (like in case of deterministic derivations of the stochastic, 
Brownian or Ornstein-Uhlenbeck type evolutions, |2J, |25|, |l7f). 

It is clear that the Madelung decomposition ip(x,t) = [exp(R + iS)](x,t) of 
a nonzero (on its domain of definition) solution of the Schrodinger equation (we 
maintain the notation D instead of h/2m, and consider the conservative case V = 
V(x)): 

id t i[>(x, t) = -DAi/>(x, t) + -L-V{x)i/)(x, t) (1) 

2mV 

whji(x,t) = DAiJi(x,t) ^— V(x)ij(x,t) 

2mD 

where tp = exp(R — iS) is a complex conjugate of ip, while R(x,t), S(x,t) are real 
functions and ip(x, 0) is taken as the initial Cauchy data for equations (1), implies 
the validity of the coupled pair of nonlinear partial differential equations 

d t p(x,t) = -V(vp)(x,t) 

Q(x,t)-V(x) = 2mD[d t S + D(VS) 2 }(x,t) (2) 

where 

p(x,t) = ip(x,t)ip(x,t) = [exp(2R)](x,t) 
v(x,t) = 2DVS(x,t) 

Q{x,t) = 2mD 2 — H —{x,t) (3) 

If, instead of complex functions we introduce real functions (we follow the 
notation of previous publications 0, H 

9(x,t) = [exp(R + S)](x,t) 

6 m (x,t) = [exp(R-S)](x,t) (4) 

then equations (2) can be replaced by the (nonlinearly coupled via Q(x,t)) pair of 
time adjoint parabolic equations for 9(x,t),9*(x,t): 

d t e. = DAe*--±-(2Q-V)9* 
2mV 
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d t 9 = -DA9 + -±-(2Q-V)9 



(5) 



with the Cauchy data 9(x, 0), 9*(x, 0) fixed by the previous Madehing exponents 
R(x, 0), S(x, 0). In turn, they imply the validity of the equations (2) with p(x,t) = 
[exp(2R))(x,t) = 9(x,t)9*(x,t). 



On the other hand, let us notice that the adjoint pair of the Schrodinger equations 
(1) comes out || ||, 0] as a direct result of an analytic continuation in time of the 



temporally adjoint parabolic problem (call it Euclidean [|T|, [T 

d t G*(x,t) = DAe*(x,t) - -L-V(x)e*(x,t) 

2mD 

d t e(x, t) = -DAQ(x, t) + -J—V(x)Q(x, t) (6) 

with a suitable (the same as in (1)) potential function V(x), defining a holomor- 
phic semigroup exp(—tH/2mD),t > 0, and thus a consistent system of solutions 
0*(x, t), Q(x,t), which gives rise to the probability measure with the quantally fac- 
torized density (BB*)(x, t), on all (finite) time intervals run by the time parameter 
t. Its relevance for the standard nonequilibrium statistical physics processes, we 
have discussed elswhere, || ||. 

An indirect effect of the analytic continuation is the mapping of the parabolic 
system (6) into rather complicated (nonlinear coupling) parabolic system (5) which 
is a mathematical, eventually probabilistic, equivalent of the Schrodinger equation. 
The seemingly strange form of (5) does not preclude the full-fledged stochastic anal- 
ysis. In fact, the standard methods appropriate for the problem (6) need only a 
slight generalization to encompass the time- dependent potentials, and next some 
boundary data analysis to deal with the (a priori admitted by (1)) nodal surfaces 
of the probability distribution, see e.g. |TT], ||, f|, [5], |37|]. Albeit, in case of (5) 
and (6), if (1) is not invoked at all, the time adjoint parabolic equation might look 
annoying for the reader unfamiliar with the properties of fundamental solutions of 
the parabolic equations (assuming their existence in the present context, cf. @, ||). 
Let us stress that there is no conflict with the traditional intuition about physically 
irreversible random transport phenomena. 



In each of the considered problems (1), (5), (6), the probability density was nat- 
urally associated with the temporally adjoint pair of partial differential equations. 
Let us choose a concrete time interval t G [0,T] and consider the boundary data 
p(x, 0), p(x, T) of the respective probability density, which we demand to be strictly 
positive on their domain of definition. We are interested in deducing a stochastic 
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process taking place in this time interval, which either induces a continuous propaga- 
tion (is measure preserving in the stationary case) of a probability density between 
the boundary data, or is consistent with the time evolution of p(x,t),t G [0,T], if 
given a priori as in case of (1) and (5). 

Since the global existence/uniqueness theorems ||, |23|, ^TJ tell us that the perti- 



nent processes might be Markovian, we fall into the well established framework, 
where for any two Borel sets A, B C R on which the respective strictly posi- 
tive boundary densities p(x, 0) and p(x, T) are defined, the transition probability 
m(A, B) from the set A to the set B in the time interval T > has a density given 
in a specific factorized form: 

m(x,y) = f(x)k(x,0,y,T)g(y) 
m(A,B)=! dxf dym( X ,y) 

J A JB 

dym(x, y) = p(x, 0) , / dxm(x, y) = p(y, T) (7) 



Here, f(x),g(y) are the a priori unknown functions, to come out as solutions 
of the integral (Schrodinger) system of equations (7), provided that in addition to 
the density boundary data we have in hands any strictly positive, continuous in 
space variables function k(x,0,y,T). Our notation makes explicit the dependence 
(in general irrelevant) on the time interval endpoints. It anticipates an important 
restriction we shall impose, that k(x, 0, y, T) is a particular form of a strongly contin- 
uous dynamical semigroup kernel: it will secure the Markov property of the sought 
for stochastic process. 

It is the major mathematical discovery 0, Q that the Schrodinger system (7) 
of integral equations admits a unique solution in terms of two nonzero, locally in- 
tegrable (i.e. integrable on compact sets) functions f(x),g(y) of the same sign 
(positive, everything is up to a multiplicative constant). 

If k(y,0,x,T) is a particular, confined to the time interval endpoints, form of 
a concrete semigroup kernel k(y, s, x, t), < s < t < T, let it be a fundamental 
solution associated with (6) (whose existence a priori is not granted), then there 
exists [§, |, |, 0, a function p(y, s, x, t): 

p(y, s, x, t) = k(y, s, x, t) ® Xj (8) 

where 

9(x,t) = dyk(x,t,y,T)g(y) (9) 
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6*(y,s)= / dxk(x,0,y,s)f(x) 



which implements a consistent propagation of the density p(x,t) = 9(x,t)8*(x,t) 
between its boundary versions, according to: 

p(x,t) = J p(y,s,x,t)p{y,s)dy (10) 

< s < t < T 

For a given semigroup which is characterized by its generator (Hamiltonian), the 
kernel k(y,s,x,t) and the emerging transition probability density p(y,s,x,t) are 
unique in view of the uniqueness of solutions f(x),g(y) of (7). For Markov processes, 
the knowledge of the transition probability density p(y, s, x, t) for all intermediate 
times < s < t < T suffices for the derivation of all other relevant characteristics. 

In the framework of the Schrodinger problem the choice of the integral kernel 
k(y,0,x,T) is arbitrary, except for the strict positivity and continuity demand. As 
long as there is no "natural" physical motivation for its concrete functional form, the 
problem is abstract and of no direct physical relevance. However, in the context of 
parabolic equations (5) and (6), this "natural" choice is automatically settled if the 
Feynman-Kac formula can be utilized to represent solutions. (Notice that in case of 
(5) the finite energy condition /[(V'?/'! 2 + V(x)p(x, t)]dx < oo, p = \?p\ 2 secures the 
boundedness from below of the potential 2Q(x, t) — V(x), see e.g. [[| [H], Effi]). 

Indeed, in this case an unambigous strictly positive semigroup kernel which is a 
continuous function of its arguments, can be introduced for a broad class of (admis- 
sible potentials. Time dependent potentials are here included as well. More- 



over, in Ref. || we have discussed a possible phenomenological significance of the 
Feynman-Kac potentials, as contrasted to the usual identification of Smoluchowski 
drifts with force fields affecting particles (up to a coefficient) in the standard theory 
of stochastic diffusion processes. 



There is an enormous literature on the Kac integral kernel issue [|Tt], [TT], [^J based 
on the concept of the conditional Wiener measure. Let us however mention that 
strictly positive semigroup kernels generated by Laplacians plus suitable potentials 
are very special examples in a surprisingly rich encompassing family. The concept 
of the "free noise", normally characterized by a Gaussian probability distribution 
appropriate to a Wiener process, can be extended to all infinitely divisible proba- 
bility distributions via, the well known to probabilists and mathematical physicists, 
Levy-Khintchine formula, see for example M. It allows to expand the framework 
from continuous diffusion processes to jump or combined diffusion-jump propaga- 
tion scenarios which are not necessarily Gaussian. All such (Levy) processes are 
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associated with the strictly positive dynamical semigroup kernels and the same per- 
tains to a number of cases when the free generator (minus Laplacian in the "normal" 
situation) acquires a potential term, to form a nontrivial Hamiltonian of a physical 
problem. 



In the existing probabilistic investigations || 0, |5], 0, based on the exploitation 
of the Schrodinger problem strategy, much stronger demand than any previous one 
was in use: guided by the observation that k(y,s,x,t) must be a function to allow 
for all advantages of (7), it was generally assumed that the kernel actually is a fun- 
damental solution of the parabolic equation. It means that the kernel is a function 
with continuous derivatives: first order-with respect to time, second order-with re- 
spect to space variables. Then, the transition probability density defined by (8) is 
a fundamental solution of the Fokker-Planck (second Kolmogorov) equation in the 
pair x, t of variables, and as such is at the same time a solution of the backward 
(first Kolmogorov) equation in the pair y, s. This feature was exploited in |3], || §• 

There is a number of mathematical subtleties involved in the fundamental so- 
lution notion, since in this case, the Feynman-Kac kernel must be a continuously 
differentiable function, and a solution of the parabolic equation itself. In fact, for 
suitable (not too bad) potentials, each fundamental solution of the parabolic equa- 



tion has the Feynman-Kac representation, |27|, and is both strictly positive and 



continuous integral kernel [10, IT]. The inverse statement is generally incorrect: 
Feynman-Kac kernels may have granted the existence status, even as continuous 
functions [|ll|, but may not be differentiable, and need not to be solutions of 



any conceivable partial differential equations. Even, if the Feynman-Kac path in- 
tegral representation applies to explicit solutions of the parabolic equations, which 
are generated from the smooth initial data by the strongly continuous semigroup 
action of the type [exp(— tH/2mD)f](x) = 9*(x,t), compare e.g. Eq.(9). 

To our knowledge, this complication in the study of Markovian representations 
of the Schrodinger interpolating dynamics (and the quantum Schrodinger picture 
dynamics in particular) for the first time has been addressed and solved in Ref . || . 

As well, the subject of the (continuous) differentiability of Feynman-Kac ker- 
nels seems to have been left aside, also in the specialized monographs [T(], [H], |27| . 
Nevertheless, we can firmly repeat the conclusion of our previous paper || that to 
give a definite (unique) Markov solution of the Schrodinger stochastic interpolation 
problem, in particular for the case of the Schrodinger picture quantum dynamics (1), 
a suitable (compatible with (5)) Feynman-Kac semigroup with its strictly positive 
and continuous in all variables kernel must be singled out. As it appears, the kernel 
may not be a fundamental solution of a parabolic equation. Anyway, for each chosen 
kernel, the associated Markov process is defined uniquely by (7)-(8), though not in 
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reverse. 



2 When diffusion processes ? 

The strategy of deducing a probabilistic solution of the Schrodinger boundary data 
problem in terms of Markov stochastic processes running in continuous time, was 
accomplished in Ref. in a number of steps accompanied by the gradual strength- 
ening of restrictions imposed on the Feynman-Kac potential. To deal with the 
commonly accepted diffusion process notion, certain additional restrictions need to 
be imposed to guarrantee that the mean and variance of the infinitesimal displace- 
ments of the continuous process have the standard meaning of the drift and diffusion 
coefficient, respectively, f32fl . 

According to the general wisdom, diffusions arise in conjunction with the parabolic 
evolution equations, since then only the conditional averages are believed to make 
sense in the local description of the dynamics. It is not accidental that forward 
parabolic equations (6) are commonly called the generalized diffusion equations. 
Also, the fact that the Feynman-Kac formula involves the integration over sam- 
ple paths of the Wiener process, seems to suggest some diffusive features of the 
Schrodinger interpolation, even if we are unable to establish this fact in a canonical 
manner. 

Clearly, the conditions valid for any e > 0: 

(a) there holds lim t i s ^ J\ y _ x \ >€ p(y, s, x, t)dx = 0, (notice that (a) is a direct con- 
sequence of the stronger, Dynkin condition, (34)), 

(b) there exists a drift function b(x, s) = lim t ± s ^ f\ y _ s \ <e (y — x)p(x, s, y, t)dy, 

(c) there exists a diffusion function a(x, s) = Um t i s -j-^ Ji y _ x i <e (y — x) 2 p(x, s, y, t)dy, 
are conventionally interpreted to define a diffusion process, ]32| . 

If we exploit the propagation formula for p(x,t), (10) and ask for the circum- 
stances under which p(x, t) is a solution of a suitable (Fokker-Planck) parabolic 
differential equation, it appears (see e.g. chap. 4.4 in Ref. [j3!| ) , that the above con- 
ditions (a), (b), (c) appear to be sufficient but not necessary to achieve this goal. 
Obviously, they can be satisfied if p(y, s, x, t) is a fundamental solution. 



Usually, one accepts that sample paths of the Wiener process are continuous 

(proposal, ac- 



21, 2 



with probability one and makes a kinematical assumption 
cording to PBfl ) by considering only these processes with continuous trajectories 
which can be derived by suitable modifications of the Wiener noise, and thus are 
regarded as being of diffusive type from the beginning. It is at this point, where the 
modern theory of stochastic differential equations and related probability measures 
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intervenes, [57], P|, p|, P^j. Then, an absolute continuity of measures relative to the 
Wiener one, allows for a continuous (and eventually diffusion process) realization 
of the Schrodinger interpolation problem. It arises in terms of weak (since an ini- 
tial probability density po(x) is attributed to the random variable X(t)) solutions 
X(t) = J" b(X(s), s)ds + v^2 W(t) of respective stochastic differential equations. 
Here, W(t) stands for the standard Wiener noise, and b(x, s) is a forward drift of 
the diffusion process. Rules of the stochastic Ito calculus allow to deduce the partial 
differential (Fokker Planck or second Kolmogorov) equation governing the dynamics 
of the probability density (and of the transition density in particular) associated 
with the process. 

Since, in the present framework, the Feynman-Kac semigroup kernel and the 
related parabolic equations (5), (6) are the principal building blocks for all the 
derivations, it seems instructive to indicate the standard procedures,^, [|, linking 
parabolic equations with diffusion processes. All of them are based on the exploita- 
tion of fundamental solutions. 

We take for granted the Feynman-Kac representation of the continuous and 
strictly positive kernel associated with the forward parabolic equation. Let us as- 
sume that we have given a bounded solution u(x, t) = J k(y, 0, x, t)u(y, 0)dz of (5) or 
(6). Let us consider u(x,t) and u(x, t + At), < At <C 1. Since u(x,t) is a solution, 
we have granted the existence of the time derivative and the validity of Taylor series 
with respect to As, at least to the second expansion order. The same (at least to 
third expansion order) applies to the Taylor expansion of u(y, t) = u(x + (y — x),t) 
about x. Consequently: 

u(x, t + At) — J k(y,t,x,t + At)u(y, s)dy ~ (11) 

k(y,t,x,t + At)[u(x,t) + (y - x)V x u(x,t) + ^(y - x) 2 A x u(x,t) + ...}dy 
On the other hand, we have 

u(x,t + At) ~u(x,t) + d t u(x,t)At (12) 

and an obvious expansion, in terms of moments of the kernel k(y, s, x, t) does emerge: 

d t u(x, t) At ~ u(x, t + At) - u(x, t) ~ -u(x, t) [1 - f k(y, t,x,t+ At)dy]+ (13) 



[V x u(x,t)\ J (y-x)k(y,s,x,t+At)dy+[—A x u(x,t)} J (y-x) 2 k(y,t,x,t + At)dy+... 
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Clearly, to reconcile this expansion with the forward parabolic equation obeyed by 
u(x,t), i.e. dtu = —cu + Am, one needs to verify whether the correct limiting 
properties are respected by the Feynman-Kac kernel. In case they would hold true, 
the arguments of Ref. would convince us that we are dealing with the diffusion 
process. 

Presently, || , the rigorous demonstration is available in case, when the kernel is 
not a fundamental solution of the parabolic equation. 

3 From positive to nonnegative solutions of para- 
bolic equations 

Following Refs. || |9|], let us focus our attention on stochastic Markov processes 
of diffusion-type (see Ref. ]7| for a jump process alternative), which are associated 
with the general temporally adjoint pair of parabolic partial differential equations: 

d t u(x, t) = Am(s, t) — c(x, t)u(x, t) (14) 

d t v(x, t) = —Av(x, t) + c(x, t)v(x, t) 

Here, c(x, t) is a real function (left unspecified at the moment) and the solutions 
u(x, t), v(x, t) are sought for in the time interval [0, T] under the boundary conditions 
set at the time-interval borders: 



Po(x) = u(x, 0)v (x, 0) (15) 
Pt(x) = u(x,T)v(x,T) 
p {x)dx = p {A) , / p T (x)dx = p T {B) 

A JB 

We assume that p is a probability measure with the density p(x), and A, B stand 
for arbitrary Borel sets in the event space. In the above, suitable units were chosen 
to eliminate inessential in the present context (dimensional) parameters, and the 
process is supposed to live in/on R 1 . 

As emphasized in the previous publications, the key ingredient of the 

formalism is to specify the function c(x, t) such that exp[— / ' H{r)dr] can be viewed 
as a strongly continuous semigroup operator with the generator H(t) = —A + c(t), 
associated with the familiar |TT| Feynman-Kac kernel: 



(f,exp[- / H{r)dr\g) = dy dxf(y)k(y,0,x,t)g(x) = (16) 
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/(w(0))£(u;(f))exp[- / c(uj(r),r)dT}dfi (uj) 
Jo 

The exponential operator should be understood as the time-ordered expression. Here 
f,g are complex functions, u(t) denotes a sample path of the conventional Wiener 
process and dpo stands for the Wiener measure. Clearly, the kernel itself can be 
explicitly written in terms of the conditional Wiener measure dpy^ s l pinned at space- 
time points (y, s) and (x, t), < s < t < T: 

k(y, s, x,t) = J exp[- j c(w(r), r)dr] dfi^(u) (17) 

As long as we do not impose any specific domain restrictions on the semigroup 
generator H(t), the whole real line R 1 is accessible to the process. Various choices 
of the Dirichlet |TT , |J boundary conditions can be accounted for by the formula (3). 



If we replace R l by any open subset Q G R 1 with the boundary dQ, it amounts to 
confining Wiener sample paths of relevance to reside in (be interior to) Q, which in 
turn needs an appropriate measure d$fl(u G ST) in (4). This is usually implemented 
by means of stopping times for the Wiener process, JT5], f|, [371 . 

Let f(x), g(x) be two real functions such that: rriT(x,y) = f(x)k(x,0,y,T)g(y) 
defines a bi-variate density of the probability measure, i.e. a transition probability 
of the propagation from the Borel set A to the Borel set B to be accomplished in 
the time interval T. In particular, we need the marginal probability densities to 
be defined: po(x) = m^(x, SI) , priy) = ^riSl^y) where Q C -R 1 is a spatial area 
confining the process. 

These formulas can be viewed as special cases of (7), so establishing an apparent 
link between the Schrodinger problem and the Feynman-Kac kernels, together with 
the related parabolic equations. Assuming that marginal probability measures and 
their densities are given a priori, and a concrete Feynman-Kac kernel (17) (with or 
without Dirichlet domain restrictions) is specified, we are within the premises of the 
Schrodinger boundary data problem. 

Let Q = Q U dQ be a closed subset of R , or R 1 itself. For all Borel sets (in 
the cr-field generated by all open subsets of TT) we assume to have known po(A) and 
Pt{B), hence the respective densities as well. If the integral kernel k(x,0,y,T) in 
the expression (5) is chosen to be continuous and strictly positive on Q, then the 
integral equations (7) can be solved Q with respect to the unknown functions f(x) 
and g(y). The solution comprises two nonzero, locally integrable functions of the 
same sign, which are unique up to a multiplicative constant. 

If, in addition, the kernel k(y, s, x, t), < s < t < T is a fundamental solution of 
the parabolic system (14) on R 1 (i.e. is a function which solves the forward equation 
in (x, t) variables, while the backward one in (y, s)), then we have defined a solution 
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of the system (14) by: 



u(x, t) = f(x, t) = / f(y)k{y, 0, x, t)dy (18) 



v(x,t) = g{x,t) = / k(x,t,y,T)g(y)dy 



Moreover, p(x,t) = f(x,t)g(x,t) is propagated by the Markovian transition proba- 
bility density: 

g (x, t) 

p(y,s,x,t) = k(y,s,x,t) ' (19) 



p(x, t) = J p(y, s)p(y, s, x, t)dy 
0<s<t<T 
d tP = Ap- V(bp) 

the result, which covers all traditional Smoluchowski diffusions ^,0. In that case, 
c(x, t) is regarded as time- independent, and the corresponding stochastic process is 
homogeneous in time. The Dirichlet boundary data can be implemented as well, thus 
leading to the Smoluchowski diffusion processes with natural boundaries, 0. Then, 
k(y, s, x, t) stands for an appropriate Green function of the parabolic boundary-data 
problem, with the property to vanish at the boundaries dQ of Q. 

Let us mention that for time-independent potentials, c(x, t) = c(x) for all t G 
[0,T], a number of generalizations is available |15|, || f|, |37|]-|39| to encompass 
the nodal sets of p{x) and hence of the associated functions f(x),g(x). The drift 
b(x) = Wlnp(x) = singularities do not prohibit the existence of a well defined 
Markov diffusion process (9), for which nodes are unattainable. In the considered 
framework they are allowed only at the boundaries of the connected spatial area Q 
confinig the process. 

The problem of relaxing the strict positivity (and/or continuity) demand for 
Feynman-Kac kernels is nontrivial |3B], with respect to the eventual construc- 



tion of the unique Markov process (9). To elucidate the nature of difficulties un- 
derlying this issue, we shall consider quantally motivated examples of the parabolic 
dynamics (14). 
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4 Nonlinear parabolic dynamics with unattain- 
able boundaries 



Let us choose the potential function c(x, t) as follows: 

x 2 1 

c(x. t) = — — (20) 

v ' ; 2(1 + t 2 ) 2 1 + t 2 K ' 

for x e R 1 , t G [0, T\. In view of its local Holder continuity (cf. Ref. ||) with expo- 
nent one, and its quadratic boundedness, the fundamental solution of the parabolic 



system is known to exist p7L HO, 41, 42]. It is constructed via the parametrix 



method. Among an infinity of regular solutions of (14) with the potential (20), we 
can in particular identify solutions of the Schrodinger boundary data problem 
for the familiar (quantal) evolution: 



,y»2 ,y» 2 

•Xj -. , , i- t O \ i 1 / 1 r t-C 



p (x) = ^y^expl--} — > p(x, t) = [2k{1 + t2 T 1/2 ™Pl-^YT¥) ] (21) 
They read 

u(x, t) = fix, t) = [2tt(1 + t 2 )r 1/4 exp(-^-^l + \ a rctan t) (22) 

4 1 + t 2 2 

x 2 1 — t 1 

v(x,t) = g(x,t) = [2tt(1 + t 2 )}~ l/4 exp(-— - -arctant) 

and, while solving the nonlinear parabolic system (14) (with c = Ap 1 / 2 / p 1//2 ), in 
addition they imply the validity of the Fokker-Planck equation: 

p{x, t) = f(x, t)g(x, t) -> d t p = Ap - V(6p) 

6x,t =2^4r = -— ? i 23 
g(x, t) 1 +t z 

Notice that p(y, s, x, t) = k(y, s, x, t) ^r'x is a fundamental solution of the first and 
second Kolmogorov (e.g. Fokker-Planck) equations in the present case. 

Let us recall that a concrete parabolic system corresponding to solutions (22) 
looks badly nonlinear. Our procedure, of first considering the linear system (but with 
the potential "belonging" to another, nonlinear one), and next identifying solutions 
of interest by means of the Schrodinger boundary data problem, allows to bypass this 
inherent difficulty. In connection with the previously mentioned quantal motivation 
of ours, let us define g = exp(R + S), f = exp(R — S) where R(x, t), S(x, t) are real 



13 



functions. We immediately realize that (5), (14) provide for a parabolic alternative 
to the familiar Schrodinger equation and its temporal adjoint: 

idti) = -Aif) (24) 

idtif) = Aip 

with the Madelung factorization xjj = exp(R + iS), ip = exp(R — iS) involving the 
previously introduced real functions R and S. 

Things seem to be fairly transparent when the parabolic system (5) or (6) allows 
for fundamental solutions. However, even in this case complications arise if nodes 
of the probability density are admitted. The subsequent discussion has a quantal 
origin, and comes from the free Schrodinger propagation with the specific choice, 
|J, of the initial data: 

2 

i> (x) = (2tt)- 1/4 x exp(~) — > (25) 

2 

Mz,t) = (2 7 r)- 1/4 (r ^ ?7 , eX ^~MT^ 

such that our nonstationary dynamics example displays a stable node at x = for 
all times . 

The parabolic system (1) in this case involves the potential function: 

V ' ; P 1/2 (x,t) 2(1 + t 2 ) 2 1 + t 2 K J 

p(x,t) = (2tt)- 1 / 2 (1 +t 2 )- 3 / 2 x 2 exp[ ' 



2(l+t 2 ) J 

The polar (Madelung) factorization of Schrodinger wave functions implies: 

R(x,t) = lnp 1/2 (x,t) (27) 

x > — > S(x, t) = S+(x, t) = — — -arctan t 

x < — > b[x, t) = b-[x,t) = — — -arctan t + n 

4 1 + t 2 2 

Although S(x,t) is not defined at x — 0, we can introduce continuous functions 
/ = exp(R — S) and g = exp(R + S) by employing the step function e(x) = if 
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x > and e(x) — 1 if x < 0. Then, the candidates for solutions of the parabolic 
system (5) with the potential (20) would read: 

v(x, t) = g(x, t) = (27r) -1 / 4 (l + t 2 )~ 3 / 4 |x| exp( — — - — ^) exp[— -arctan t + 7re(x)] 

+ 1 ^ (28) 

x 2 1 + t 3 

■u(x, i) = f(x, t) = (27r) _1//4 (l + t 2 )~ 3 ^ 4 \x\ exp(— — -) exp[-arctan t — ire(x)] 

4 1 ~\~ t 2 

For all x ^ we can define the forward drift 

b(x,t) =2 , v '/ = - -x- r 29 

which displays a singularity at x — 0. Nonetheless, (bp)(x,t) is a smooth function 
and the Fokker-Planck equation <9 t p = Ap — V(fep) holds true on the whole real 
line R 1 , for all t e [0,T]. Notice that there is no current through x — 0, since 
v(x,t) = 2VS(x,t) = vanishes at this point for all times. 

Our functions f(x, t),g(x, t) are continuous on R 1 , which however does not imply 
their differentiability. Indeed, they solve the parabolic system (5) with the potential 
(20) not on R 1 but on (— oo,0) U (0, +oo). Hence, almost everywhere on R 1 , with 
the exception of x = 0. 

An apparent obstacle arises because of this subtlety: these functions are not even 
weak solutions of (1), because of: 

/+oo r+oo 
d t f(x,t)(f)(x)dx+ / Vf(x,t)V(f)(x)dx+ (30) 

1 r+oo 

- / c(x,t)f(x,t)<j)(x)dx ^ 

2 J —oo 

for every test function such that 0(0) 7^ 0, continuous and with support on a 
chosen compact set (e.g. vanishing beyond this set). 

One more obstacle arises, if we notice that c(x,t), (20) permits the existence 
of the unique, bounded and strictly positive fundamental solution for the parabolic 
system (14). Then, while having singled out a fundamental solution and the bound- 
ary density data po(x), Pt{%) consistent with (20), we can address the Schrodinger 
boundary data problem associated with (2), (3): 

u(x, 0) J k(x, 0,y,T)v(y,T)dy = p (x) (31) 
v(x,T) J k(y,0,x,T)u{y,0)dy = p T {x) 
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expecting that a unique solution u(x, 0),v(x,T) of this system of equations implies 
an identification u(x, 0) = f(x, 0) and v(x, T) = g(x, T). 

However, it is not the case and our f(x,t),g(x,t) do not come out as solutions 
of the Schrodinger problem, if considered on the whole real line R 1 , on which the 
fundamental solution sets rules of the game. Indeed, let us assume that (31) does 
hold true if we choose u(x, 0) = f(x, 0), v (x, T) = g(x,T), with / and g defined by 
(21). Since, in particular we have 

g(x, T) J k(y, 0, x, T)f(y, 0)dy = g(x, T)f(x, T) (32) 

then for there holds: 

f(x,T) = J k(y,0 } x,T)f(y,0)dy (33) 

Both sides of the last identity represent continuous functions, hence the equality is 
valid point-wise (i.e. for every x). We know that f(y, 0) is continuous and bounded 
on R 1 , and k(y, 0,x, T) is a fundamental solution of (1). Hence the right-hand-side 
of (33) represents a regular solution of the parabolic equation. Such solutions have 
continuous derivatives, while our left-hand-side function f(x,T) certainly does not 
share this property. Consequently, our assumption leads to a contradiction and (33) 
is invalid in our case. 

It means that the fundamental solution (e.g. the corresponding Feynman-Kac 
kernel) associated with (16) is inappropriate for the Schrodinger problem analysis, 
if the interpolating probability density is to have nodes (i.e. vanish at some points). 

In our case, x = is a stable node of p(x,t), and is a time-independent re- 
pulsive obstacle for the stochastic process. An apparent way out of the situation 
comes by considering two non-communicating processes, which are separated by the 
unattainable barrier at x = 0, ]57], f£3"| . The pertinent discussion can be found 

in 1. 

5 The "Wiener exclusion" 

The conventional definition of the Feynman-Kac kernel (in the conservative case) 

exp[-t(-A + c)](y,x) = J exp[- j\{u{r))dr]d^{uj) (34) 

comprises all sample paths of the Wiener process on R l , providing merely for their 
nontrivial redistribution by means of the Feynman-Kac weight exp[— Jq c{uj{r))dr] 
assigned to each sample path uj{s) : u>(0) = y,uj(t) = x. 
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Assume that c(x) is bounded from below and locally (i.e. on compact sets) 
bounded from above. Then, the kernel is strictly positive and continuous |TT| . 
For c = we deal with the conditional Wiener measure 

exp(tA)(y,x) = /$$Ks) G R 1 ; < s < t] = (35) 

fi[u(s) G R 1 ; u(0) = y, u(t) = x;0<s<t] 

pinned at space-time points (y, 0) and (x, t). 

The previous discussion indicates that _R_ is inaccessible for all sample paths 
originating from R + . In reverse, R + is inaccessible for those from R_. As well, we 
may confine the process to an arbitrary closed subset Q C R 1 , or enforce it to avoid 
("Wiener exclusion" of Ref. |[30|| ) certain areas in R 1 . 

In this context, it is instructive to know that [jnj for an arbitrary open set Q, 
there holds: 

expit A n ) (y, x) = ^ [u(s) G Q, < s < t) (36) 

which is at the same time a definition of the operator — Aq, i.e. the Laplacian with 
Dirichlet boundary conditions, and that of the associated semigroup kernel. This 
formula provides us with the conditional Wiener measure which is confined to the 
interior of a given open set, ^3, ^5], [5|. 



We can introduce an analogous measure, which is confined to the exterior of a 
given closed subset S C R 1 . In case of not too bad sets (like an exterior of an 
interval in R l or a ball in R n , the corresponding integral kernel in known |L1] to 



be positive and continuous. Technically, if S is a (regular) closed set such that the 
Lebesgue measure of dS is zero, then: 

exp(tA R \ S )(y, x) = ^%[oo{s) ^S;0<s<t] (37) 

The Feynman-Kac spatial redistribution of Brownian paths can be extended to 
cases (36), (37) through the general formula valid for any f,g G L 2 (fl), where f2 is 
any open set of interest (hence R\S, in particular): 

(f,exp(-tH n )g) = [ 7{u{ti))g{u{t))exp[- f c{uj{r))dr}d^{uj) (38) 
Jci Jo 

It gives rise to the integral kernel comprising the restricted Wiener path integration, 
which is defined at least almost everywhere in x, y. Then, its continuity is not 
automatically granted. We can also utilize a concept of the first exit time Tq for the 
sample path started inside fl (or outside S) 

T n (u) = inf[t > 0,X t (u) i n] (39) 
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where X t is the random variable of the process. Then, we can write, 0, |5[] 

exp{-tH n ){y, x) = J exp[- j c(u(r))dT}d/i^[u; t < T n ] = (40) 

J n ex P{~ J Q c M T )) dr ]^(S)V) 

It is an integration restricted to these Brownian paths, which while originating from 
y G Q at time t = are conditioned to reach x G Q at time £ > without crossing 
(but possibly touching) the boundary dS of S. The contribution from paths which 
would touch the boundary without crossing, for at least one instant s G [0,t] is of 
Wiener measure zero, 

In case of processes with unattainable boundaries, with probability 1, there is 
no sample path which could possibly reach the barrier at any instant s < oo. 

The above discussion made an implicit use of the integrability property 

/ c(u(s))ds < oo (41) 
Jo 

for uj G R l , < s < t., in which case the corresponding integral kernel (for bounded 
from below potentials) is strictly positive. Then, if certain areas are inaccessible to 
the process, it occurs excusively |22|]-|39| due to the drift singularities, which are 
capable of "pushing" the sample paths away from the barriers. 

The previous procedure can be extended to the singular [fj6|-|54l potentials, 
which are allowed to diverge. Their study was in part motivated by the so called 
Klauder's phenomenon (and the related issue of the ground state degeneracy of 
quantal Hamiltonians), and had received a considerable attention in the literature. 

In principle, if S is a closed set in R 1 like before, and c(x) < oo for all x G 
Q = R\S, while c(x) = oo for x G S, then depending on how severe the singularity 
is, we can formulate a criterion to grant the exclusion of certain sample paths of 
the process and hence to limit an availability of certain spatial areas to the random 
motion. Namely, in case of (41) nothing specific happens, but if we have 

/* c(u(r))dT = oo (42) 
Jo 

for oj(t) G S for some r G [0,t], then the "Wiener exclusion" certainly appears: we 
are left with contributions from these sample paths only for which (42) does not 
occur. Unless the respective set is of Wiener measure zero. 

The area Q comprising the relevant sample paths is then selected as follows: 

Q, = b; f ' c(u(r))dT < oo] (43) 
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In particular, the criterion (43) excludes from considerations sample paths which 
cross S and so would establish a communication between the distinct connected 
components of Q. 

The singular set S can be chosen to be of Lebesgue measure zero and contain a 
finite set of points dividing R into a finite number of open connected components. 
With each open and connected subset Q C -R 1 we can |46| associate a strictly positive 
Feynman-Kac kernel, which can be expected to display continuity. 

Since the respective potentials diverge on S, their behaviour in a close neigh- 
bourhood of nodes is quite indicative. For, if u>s is a Wiener process sample path 
which is bound to cross a node at < s < t, then the corresponding contribution to 
the path integral vanishes. Such paths are thus excluded from consideration. If their 
subset is sizable (of nonzero Wiener measure), then the eliminated contribution 

/ exp[- f c(u s (r))dr]dfi^(u) = (44) 

is substantial in the general formula (35). 

At the same time, we get involved a nontrivial domain property of the semigroup 
generator H = — A + c resulting in the so called ground state degeneracy j|B|, E7J f4"g| . 
Let us recall (Theorem 25.15 in Ref. |l 1]]) that if c is bounded from below and locally 
bounded from above, then the ground state function of H = —A + c is everywhere 
strictly positive and thus bounded away from zero on every compact set. 



6 Singular potentials, ground state degeneracy and 
the "Wiener exclusion" 

Our further discussion will concentrate mainly on singular perturbations of the har- 
monic potential. Therefore, some basic features of the respective parabolic problem 
are worth invoking. The eigenvalue problem (the temporally adjoint parabolic sys- 
tem now trivializes): 

-Ag + (x 2 — E)g = = Af — (x 2 - E)f (45) 

has well known solutions labeled by E n = 2n + 1 with n — 0, 1, 2, .... In particular, 
go(x) = fa(x) = -^j^expi—^) is the unique nondegenerate ground state solution. 
The corresponding Feynman-Kac kernel reads 

exp(-tH)(y,x) = k(y,0,x,t) = k t (y,x) = 
(vr)-^! _ ex P (-t))^ 2 ex P [-^^ - (y^Pi-t)-*) 2 } (46 ) 
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d t k = -A x k + {x 2 - l)k 

and the invariant probability density p(x) = f(x)g(x) = (n)~ 1 / 2 exp(—x 2 ) is pre- 
served in the course of the time-homogeneous diffusion process with the transition 
probability density 

g (x) 

p(y, s, x, t) = k t _ s (y, x)—— (47) 

9W 

We have p(y, s, x, t) = p(y, 0, x, t — s). 

Notice the necessity of the eigenvalue correction (renormalization) of the poten- 
tial, both in (34) and (38), which is indispensable to reconcile the functional form 
of the forward drift b(x) = 2Vlng(x) = —2x with the general expression for the 
corresponding (to the diffusion process) parabolic system potential 

1 b 2 

c = c(x,t)=dtlng + -(- + Vb) (48) 

which equals c(x) = x 2 — 1 in our case. 

Let us pass to the singular (degenerate) problems. 
The canonical (in the context of Refs. |46[]-[^l|]) choice of the centrifugal potential: 

c E (x) = x 2 + ^-E (49) 



x 



generates a well known spectral solution 0, ^TJ for Hg = [—A + CE(x)]g. The 
eigenvalues: 

£ n = 4n + 2 + (l + 8 7 ) 1/2 (50) 

with n = 0,1,2,... and 7 > — ~ =^ (1 + 87) 1 / 2 > y/2, are associated with the 
eigenfunctions of the form: 

g n {x)=x^' 2 exp{-^)L«{x 2 ) (51) 



a = (1 + 87) 1/2 

(n + a)\ {-x 2 f 
u=0 (n - u)\(a + v)\ v\ 

L%{x 2 ) = 1 , L^(x 2 ) = -x 2 + a + l 

It demonstrates an apparent double degeneracy of both the ground state and of the 
whole eigenspace of the generator H. The singularity at x = does not prevent 
the definition of H = —A + x 2 + ^ since this operator is densely defined on an 
appropriate subspace of L 2 (R l ). This singularity is sufficiently severe to decouple 
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(— oo, 0) from (0, oo) so that L 2 (— oo, 0) and L 2 (0, oo) are the invariant subspaces of 
H with the resulting overall double degeneracy. 



Potentials of the form, |3~0|, pfl: 



c( x ) = x 2 + [dist(x, dtt)}- 3 (52) 

where dfl can be identified with dS, and S is a closed subset in R 1 of any (zero or 
nonzero) Lebesgue measure, have properties generic to the Klauder's phenomenon. 
Because the Wiener paths are known to be Holder continuous of any order ^ — e, e > 
and of order | in particular, there holds J * c(u(t))cIt = oo if uj{r) E S for some r. 
Conversely, J c(co(r))dT < oo if to never hits S. This implies that the relevant 
contributions to: 

(/, exp[-t(-A + c)\g) = J J(u(0))g(u(t))exp[- jf* c{u{r))dr\d^{u) (53) 

come only from the subset of paths defined by 

Qt = [u; f c(uj(T))dr < oo] (54) 
Jo 

The above argument might seem inapplicable to the centrifugal problem. However 
it is not so. In the discussion of the divergence of certain integrals of the Wiener 
process, in the context of Klauder's phenomenon, it has been proven |53) that for 
almost every path from x — 1 to x = — 1 (crossing the singularity point x = 0) there 
holds |a;(r)| _1 (ir = oo for any 5 > 0. 

To be more explicit: if t\ = T\(oS) is the first time such that the Wiener process 
W(t) = W(t,u) attains the level (location on R 1 ) W{t\) = 1, then the integral over 
any right-hand-side neighbourhood (n, T\ + 5) of T\ diverges: 

Tl+8 

c(u(t) - l)dt = oo (55) 

Tl 

if J^i c(x)dx = oo. 

In case of the left-hand-neighbourhood of Ti, we have 

c{u(t) - l)dt = oo (56) 

n-<5 

if xc(x)dx = oo. 

All that holds true in case of the centrifugal potential, thus proving that the 
only subset of sample paths, which matters in (54) is (55). Obviously, Q t does not 
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include neither paths crossing x = nor those which might hit (touch) x = at 
any instant. The singularity is sufficiently severe to create an unattainable repulsive 
boundary for all possible processes, which we can associate with the spectral solution 
(50),(51). 

After the previous analysis one might be left with an impression that the ap- 
pearence of the stable barrier at x = persisting for all t G [0,T], is a consequence 
of the initial data choice ipo(0) = for the involved quantum Schrodinger picture 
dynamics. In general it is not so. For example, ipo{x) — x 2 exp(—x 2 /4) which van- 
ishes at x = 0, does not vanish anymore for times t > of the free evolution. On 
the other hand, somewhat surprisingly from the parabolic (intuition) viewpoint, the 
node can be dynamically developed from the nonvanishing initial data and lead to 
the nonvanishing terminal data. 

Let us consider, [[J, a complex function: 

«,,.) = (1 + itr^eM-^} + jfjj] (57) 

2 2 

which solves the free Schrodinger equation with the initial data i/j(x, 0) = ^exp{— ^-). 
It vanishes at x = exclusively at the initial instant t = of the evolution. 
Obviously, there is nothing to prevent us from considering 

ift(x,t)=il)(x,t-a) (58) 

for a > 0. It solves the same free equation, but with nonvanishing initial data. How- 
ever, the node is developed in the course of this evolution at time t = a and instanta- 
neously desintegrated for times t > a. Here, the Schrodinger boundary data problem 
would obviously involve two strictly positive probability densities po(x) = \^(x, 0)| 2 
and Pt(x) = \^(x,T)\ 2 , T > a. It would suggest to utilize the theory M, based on 
strictly positive Feynman-Kac kernels, to analyze the corresponding interpolating 
process. However, this tool is certainly inappropriate and cannot reproduce the a 
priori known dynamics, with the node arising at the intermediate time instant. 

To handle the issue by means of a parabolic system, which we can always asso- 
ciate with a quantum Schrodinger picture dynamics, let us evaluate the potential 
c(x,t) appropriate for (5). 

In view of 

2 4 

p(x, t) = const (1 + t 2 y 5 / 2 exp[- * ] - x 2 t 2 + t 2 (l + t 2 )] (59) 

we have (while setting w 1/2 (x, t) — — x 2 t 2 + t 2 (l + t 2 )}): 

Ap 1 / 2 (x,t) 1. x „ . 9 \, 1 . 7 i , „. 

c(x > t] = = - 4 (-TT¥ + Vw) + 2<— + Aln "> = (60) 
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1 x 2 1 3x 2 -2t 2 x 1 

4(l+t 2 ) 2 ~ 2f- t 2 X 2 +t 2 {l + t 2 ) ~ 2(1 + t 2 ) + 

1 3x 2 -2t 2 I, x 3 - 2t 2 x 



2 f - t 2 X 2 + t 2 (l + t 2 ) 4 V ^ - £2 X 2 + f 2(! + t 2) ' 

The expression looks desparately discouraging, but its t J, (i.e. the initial data ) 
limit is quite familiar and displays a centrifugal singularity at x = 0: 

Ap^foO) x 2 2 5 

cM = 7« = I + ^"2 (61) 

Since the original, dimensional expression for the centrifugal eigenvalue problem 
is |3: 



; A + ^ + i)9 = Eg (62) 



S n = m[2n + l + i(l + 8 7 ) 1 / 2 ] 

with n = 0,1,..., an obvious adjustment of constants m = 1/2,7 = 1 allows to 
identify E = 5/2 as the n = eigenvalue of the centrifugal Hamiltonian H = 

"A + ^ + J, 

A peculiarity of the considered example is that it enables us to achieve an ex- 
plicit insight into an emergence of the centrifugal singularity and its subsequent 
destruction (decay) for times t > a, due to the free quantum evolution. 

2 2 

In view of the degeneracy of the ground-state eigenfunction \exp{—^-) of the 
centrifugal Hamiltonian, we deal here with the gradually decreasing communication 
between R + and _R_, which results in the emergence of the completely separated 
(disjoint) sets (— oo, 0) and (0, +oo) at t = a, followed by the gradual increase of 
the communiaction for times t > a. By "communication" we understand that the 
set of sample paths crossing x = forms a subset of nonzero Wiener measure. 

It also involves a generalisation (cf. also Refs. §|, HI HH) t° time- dependent 
Feynman-Kac kernels: 

(/, exp[- H(r)dr}g) = J f(u(s))g(u(t))exp[- jT c(uj(t), r)dr] d^{oo) (63) 

Qs,t=[w, [ c(u (t), t) dr < oo] 



The finiteness condition f 8 c(uj(t), r)dr < oo, surely does not hold true, if 5 > 
is sufficiently small, cf. (60), (61). 



23 



Let us mention that some interesting mathematical questions were left aside in 
the present discussion. For example, even in case of conventional Feynman-Kac ker- 
nels, the weakest possible criterions allowing for their continuity in spatial variables 
are not yet established. An issue of the continuity of the kernel in case of general 
singular potentials, needs an investigation as well. The phenomenological recipes 
justifying the concrete choice of the kernel are not unequivocally established as yet, 
see e.g. our discussion of the Smoluchowski drift versus Feynman-Kac potential 
issue in Ref. [|J. 
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